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ABSTRACT 

We study magneto-elastic oscillations of highly magnetized neutron stars (magnetars) which 
have been proposed as an explanation for the quasi-periodic oscillations (QPOs) appearing in 
the decaying tail of the giant flares of soft ganmia-ray repeaters (SGRs). We extend previous 
studies by investigating various magnetic field configurations, computing the Alfven spectrum 
in each case and performing magneto-elastic simulations for a selected number of models. By 
identifying the observed frequencies of 28 Hz (SGR 1900+14) and 30 Hz (SGR 1806-20) with 
the fundamental Alfven QPOs, we estimate the required surface magnetic field strength. For 
the magnetic field configurations investigated (dipole-like poloidal, mixed toroidal-poloidal 
with a dipole-like poloidal component and a toroidal field confined to the region of field lines 
closing inside the star, and for poloidal fields with an additional quadrupole-like component) 
the estimated dipole spin-down magnetic fields are between 8 x lO^^G and 4 x lO^^G, in broad 
agreement with spin-down estimates for the SGR sources producing giant flares. A number of 
these models exhibit a rich Alfven continuum revealing new turning points which can produce 
QPOs. This allows one to explain most of the observed QPO frequencies as associated with 
magneto-elastic QPOs. In particular, we construct a possible configuration with two turning 
points in the spectrum which can explain all observed QPOs of SGR 1900+14. Finally, we 
find that magnetic field configurations which are entirely confined in the crust (if the core is 
assumed to be a type I superconductor) are not favoured, due to difficulties in explaining the 
lowest observed QPO frequencies (/ < 30 Hz). 



1 INTRODUCTION 

The discovery of quasi-periodic oscillations (QPOs) in the decay- 
ing tail of a giant flare of a soft-gamma ray repeater^ (SGR) by 
Israel et al. (2005) (see also Watts & Strohmayer 2007, and ref- 
erences therein) has stimulated strong interest in the theoretical 
modelling of oscillations of magnetized neutron stars (magnetars). 
The observed frequencies of the QPOs are roughly 18, 26, 30, 
92, 150, 625, and 1840 Hz for the outburst of SGR 1806-20 and 
28, 53, 84, and 155 Hz for SGR 1900+14, respectively Since the 
first attempts to explain these QPOs in terms of torsional, purely 
shear modes of the solid crust of the neutron star (Duncan 1998; 
Strohmayer & Watts 2005; Piro 2005; Sotani et al. 2007; Samuels- 
son & Andersson 2007; Steiner & Watts 2009) the theoretical un- 
derstanding of the oscillations has evolved. The most recent stud- 
ies (Gabler et al. 2011, 2012; Colaiuda & Kokkotas 2011b; van 
Hoven & Levin 201 1, 2012) favour global magneto-elastic oscilla- 
tions because, due to the coupling through the magnetic field, the 
core absorbs the shear modes of the crust very efficiently even for 
magnetic field strengths well below those estimated for magnetars, 

^ SGRs are assumed to be highly magnetized neutron stars showing re- 
peated bursts in the soft-gamma ray spectrum (Duncan & Thompson 1992) 



Bi5 - 1 with Bi5 = B[10^^G] (Duncan & Thompson 1992). 
Nonetheless, long-lived QPOs which have predominantly Alfven 
character have also been studied (Sotani et al. 2008; Cerda-Duran 
et al. 2009; Colaiuda et al. 2009). The frequencies of successive 
overtones of these QPOs are found in integer relations and, thus, 
can potentially explain some of the observed frequency relations. 
In two previous papers (Gabler et al. 2011, 2012) we have shown 
that these oscillations can reach the surface of the star with signifi- 
cant amplitudes only for magnetic fields stronger than Bi^ > 1, if 
the internal field has a global, dipole-like structure. 

A successful interpretation of the observed QPOs has the po- 
tential to constrain the equation of state (EoS) of the cold matter 
at supranuclear densities which occur only in the core of a neutron 
star (Samuels son & Andersson 2007). Therefore, neutron stars are 
a unique laboratory that can be used to increase our knowledge of 
the fundamental physics describing the interaction of nucleons and 
other elementary particles at the conditions inside those stars. 

While the current models of magneto-elastic oscillations are 
promising to explain some of the observed QPO frequencies, all 
existing studies are restricted to a very limited set of magnetic field 
configurations. The aim of the present work is to explore the in- 
fluence of the magnetic field configuration in the magneto-elastic 
model of QPOs to assess its validity. It must be stressed that very 
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little is known from observations about the internal magnetic field 
configurations of magnetars. The only observational constraint is 
the presence of a strong external field which is responsible for 
the observed spin-down of the magnetar. Both purely toroidal and 
purely poloidal magnetic fields are known to be unstable in un- 
stratified stars, i.e. stars with a barotropic EoS which depends on 
one parameter only (Tayler 1973; Markey & Tayler 1973). Non- 
linear simulations confirm these instabilities (Braithwaite & Spruit 
2006; Kiuchi et al. 2011; Ciolfi et al. 2011; Lasky et al. 2011; 
Lander & Jones 2011a,b) and suggest that a twisted torus config- 
uration, with mixed poloidal and toroidal field, is expected. There 
have been some attempts to model magnetic field configurations as 
equilibrium axisymmetric configurations with mixed poloidal and 
toroidal field (Colaiuda et al. 2008; Kiuchi & Kotake 2008; Ciolfi 
et al. 2009; Lander & Jones 2009). However, no stable configura- 
tion has yet been found for barotropic stars (Lander & Jones 2012). 
Nevertheless, different possibilities could stabilize magnetic fields 
in stars: a more complicated magnetic field structure (possibly non- 
axisymmetric), stratification, or the presence of a solid crust. Un- 
til this issue is settled, it is legitimate to explore the most general 
setup of magnetic field configurations possible and analyze their 
influence on the QPO frequencies, leaving aside the stability of the 
configurations. This paper, in which we explore a wide range of 
plausible axisymmetric magnetic field configurations, summarizes 
our findings in this direction. 

We begin this paper in Section 2 with a short overview of our 
theoretical and numerical framework, described in more detail in 
Cerda-Duran et al. (2009) and Gabler et al. (2012). We also de- 
scribe in detail in this section how to obtain different magnetic field 
configurations with a new implementation of the MAGSTAR tool of 
the LORENE library which we call MAGNETSTAR. In the follow- 
ing sections we consider different magnetic field configurations, in- 
cluding dipole-like fields in Section 3 and quadrupole-like fields in 
Section 4, where mixed quadrupole-dipole-like fields are also con- 
sidered. We further divide the dipole-like fields into purely poloidal 
fields (Section 3.1) and mixed toroidal-poloidal fields (Section 3.2). 
All magnetic field configurations considered up to Section 4 pene- 
trate the whole volume of the star. In Section 5 we study the case of 
a magnetic field confined to the crustal region, a situation motivated 
by the indications that neutron stars may contain superconducting 
protons (Page et al. 2011; Shternin et al. 2011) in which case the 
magnetic field could be expelled from the core of the neutron star, if 
it is a type I superconductor. A summary of our results is presented 
in Section 6. 

We use units where c — G—1 with c and G being the speed 
of light and the gravitational constant, respectively. Latin (Greek) 
indices run from 1 to 3 (0 to 3). Partial and covariant derivatives are 
indicated by a comma and a semicolon, respectively. We apply the 
Einstein summation convention. 



2 THEORETICAL FRAMEWORK 

2.1 General-relativistic magneto-hydrodynamics of elastic 
bodies 

In the present work we adopt the framework of 3+1 split of general 
relativity for a spherically-symmetric spacetime. The metric takes 
the following form 

ds^ = —a^dt^ + ^'^'^ijdx^dx'^ , 



where a is the lapse function, 
$ is the conformal factor. 



(1) 

is the flat, spatial 3 -metric, and 



The stress-energy tensor T^^ for a magnetized perfect fluid 
with elastic properties can be written as a sum of different contri- 
butions: 

^ — -'- fluid -'- magn "T J- elas 

= phu^u'' + Pg^"" + h'u^u'' + ^6'^^" - 6^6" 

-2//sE^" , (3) 

where p is the rest-mass density, h — 1 + e + p/pthe specific en- 
thalpy, e the specific internal energy, P the isotropic fluid pressure, 
the 4-velocity of the fluid and g^^ the metric tensor. In the mag- 
netic contribution to the stress-energy tensor we find the magnetic 
field 6^, measured by a co-moving observer, with 6^ = h^h^, and 
in the elastic contribution we find the shear modulus p^ and the 
shear tensor S^"^. 

We are interested in torsional oscillations of neutron stars and 
are thus allowed to simplify the model as follows: (i) We use a 
barotropic EoS. (ii) For poloidal background fields and in axisym- 
metry the axial oscillations decouple at the linear level from the 
polar oscillations. With these assumptions, we evolve only B"^ and 
S^p (see definitions below), (iii) We apply the Cowling approxima- 
tion, regarding the metric as fixed, (iv) For torsional oscillations the 
shear tensor E*-^ can be expressed as 












(4) 



while the other components vanish, = YP^ = 0. Here, we 
have introduced the displacement C » which is defined via its time 
derivative as = av^. 

With these simplifications the conservation of energy and mo- 
mentum VvT^^ — 0, baryon number conservation Vu{pu^) — 0, 
and the Maxwell's equations T^^^ = and V^F^"" = AttJ^ 
lead to the following flux-conservative hyperbolic system of equa- 
tions 



1 fd^u d^¥' 



-9 V 



: 0. 



(5) 



The various quantities introduced here are the Faraday electromag- 
netic tensor field F^'' and its dual the electric current J^, 
— OL-sJ^, and 7 = det{'jij). The state vector U and the flux 
vector F* are given by 



U = 

= 

= 



W 



(6) 
(7) 

(8) 



where W — au^ is the Lorentz factor, is the 3-velocity of the 
fluid, and Si = {ph + b^)W'^Vi — abilP is a generalization of the 
momentum density. The relation between the magnetic field mea- 
sured by an Eulerian observer, 5*, and 6* is given by 



WB'vi B' -\-W'^v^Bjv' 



W 



(9) 



Eq. (5) is complemented by the evolution equations of the dis- 
placement 



(e).*-K«),r = 0, 



(10) 

(11) 
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Boundary conditions 

The boundary condition at the neutron star surface consist of the 
continuous traction condition and the assumption that there are no 
surface currents. Therefore, the tangential magnetic field compo 
nents have to be continuous, i.e. ^^rust 

and ^^rust,r = 

At the crust-core interface we have two separate cases: (i) If 
the magnetic field penetrates the core we can define a displacement 
at this interface as well and has to be continuous. As above, 
we use the continuity of the traction to set the remaining condi- 
tion which then reads ^tore,r = [l + fis/^\b^f] ^^rust,r. (ii) If 
the magnetic field is confined to the crust, the continuous traction 
condition leads to ^ crust, r = at the crust-core interface, too. 



Numerical methods 

The numerical simulations are performed with the non-linear 
GRMHD code MCOCOA described in previous works (Cerda- 
Duran et al. 2008, 2009; Gabler et al. 2012). This code is able to 
compute the torsional oscillations of neutron stars with strong mag- 
netic fields and an elastic crust. To solve the GRMHD equations we 
employed a high-resolution shock-capturing scheme with a method 
of lines for the time-advance. 

Complementary to the numerical simulations, we will also 
make use of the semi-analytic model presented in Cerda-Duran 
et al. (2009), which allows us to compute the spectrum of standing- 
wave magneto-elastic oscillations along individual magnetic field 
lines. The model is based on an integration of a perturbation along 
the magnetic field lines in the limit of short wavelengths. We ex- 
tended this method in Gabler et al. (2012) to approximately take 
into account the elastic crust. 



Microphysics 

We construct neutron star models choosing different combinations 
of EoS: the EoS APR (Akmal et al. 1998) or the stiffer EoS L 
(Pandharipande & Smith 1975) for the core, and EoS NV (Negele 
& Vautherin 1973) or EoS DH (Douchin & Haensel 2001) for the 
crust. 

The calculation of the shear modulus /xs is based on the zero- 
temperature limit of Strohmayer et al. (1991) 



Ij^s = 0.1194 



ni{Ze) 



(12) 



where rii is the ion density, Ze the ion charge, and = 3/(47rni) 
the average ion spacing. For details see Gabler et al. (2012). 



2.2 Magnetic equilibria 

The hydrostatic equilibrium equations for a self-gravitating mag- 
netized star (see Gourgoulhon 2012, and references therein) can 
be simplified under the following assumptions: (i) axisymmetry, 
(ii) no rotation or meridional flows, (iii) no net charges, (iv) a 
barotropic EoS of the form P = P{p), and (v) an unstressed crust, 
— Q These conditions are appropriate to describe magnetars, 
which are regarded to be cold and slowly rotating. In that case 



(log h + log a) , 



(13) 



where J'^ is the magnetic current. For the magnetic field strength 
relevant for magnetars, B < 10^^ G, deformations of the equilib- 
rium configurations due to magnetic field stresses are small (Boc- 
quet et al. 1995). Therefore, we can expand around the unmagne- 
tized spherical equilibrium values, h = ho -\- h' and a = ao -\- a\ 
such that h' <^ ho, a' ^ ao, and 



(log/io + logao) i = 0. 



(14) 



Neglecting magnetic terms in density and pressure terms the result- 
ing equation for the linearized perturbations is 

1 



(15) 



Poho 

where M = h' /ho a /ao <^ 1, which quantifies the magnitude 
of the perturbation and can be used to assess the quahty of the ap- 
proximation. The equilibrium models presented later in this work 
do never exceed a maximum mas nitude of M < 5 x 10"^. The 
current and the magnetic field are related by Ampere's law 

e'^''{aBk)j=47raJ\ (16) 

Considering the restriction of the solenoidal condition for the mag- 
netic field there are only two free functions describing the mag- 
netic field configuration. For convenience we choose them to be 
Hep = aBcp/4:7T and A^p, the latter being the (y9-component of the 
vector potential Ai defined as 



The equilibrium equations as a function of A^ and read 



Mr = 



1 



Poho 
1 

poho [' 

Acp^r-Bcp^Q 



(17) 

(18) 

(19) 
(20) 



2.2.7 General magnetic field configurations 

Eqs. (18-20) lead to four different branches of possible solutions 
for the current distribution: 

• Type I.- Mixed poloidal-toroidal current (H^p = H^{A^) 
and dHcp/dA^ 7^ Oj generating a mixed poloidal-toroidal field: 



Eqs. (18-20) imply M 
1 



M{A^) and 



r = 



ip,e 



-Hin 



J"^ = poho 



-9 
dM 

dA^ 



47r7' 



ipip ^ 



dH^ 
""dA^- 



(21) 



In this case we have two free functions H^{A^) and M(A^) 
to determine the current J"^{A^), which in turn gives via the (f- 
component of Ampere's law A^. The solution of the system of 
equations is identical to that described in Bocquet et al. (1995). 
It generates a magnetic field with mixed poloidal and toroidal com- 
ponent. In the case M = 0, the configuration is force-free. 

• Type II.- Purely toroidal current (Hep = 0, A^p ^ and 
0) generating a purely poloidal field: Eqs. (18-20) imply 
M = M{A^) and 



J' = (0,0, poho 



dM 
dA^ 



(22) 
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For a given form of M{A^) and using J"^ in the v?-component of 
Ampere's law, Eq. (16), an elliptic system results that can be solved 
for A^. This is the procedure described by Bocquet et al. (1995) 
to compute self-consistent magnetic field equilibria, although their 
approach is more general since rotation, charges and the deforma- 
tions of the star are considered, too. This configuration corresponds 
to a limiting case of the type I solution for 0. 

• Type III.- External currents (H^ = 0, = and ^ 
Oj generating di force-free poloidal field: Eqs. (18-20) imply M = 
0, i.e. force-free poloidal configurations. It describes the magnetic 
field in regions with vanishing currents. The currents generating the 
magnetic field are external to the region described by this currents. 
The solution depends only on the potential A^ that is a solution of: 



V7 



0. 



(23) 



This configuration corresponds to a limiting case of a type II current 
configuration with — 0. 

• Type TV. - Purely poloidal current (A^ = and 7^ Oj gen- 
erating a purely toroidal field: Eqs. (18-20) imply M = M{H^) 
and 

dM _ 47r7'^'^ ~ 



pohoa 



(24) 



For a given M{H^), Eq. (24) provides the current and magnetic 
field configuration 



-Hin 



1 



- Hip , r 5 



(25) 



This current distribution generates a purely toroidal field confined 
to the region where J'* ^ 0. 

Once the current is known, and hence A^p and H^p, the mag- 
netic field can be computed as 



1 



1 47r-7'^'^ - 

.——A ' H 

V7 a " 



(26) 



2.2.2 Magnetic field in magnetars 

In order to describe magnetars, we make the following assump- 
tions: (i) the dipole-like component dominates the field at long 
distances and matches the observed spin down within uncertain- 
ties. This constraint excludes purely toroidal magnetic fields, (ii) 
There are no currents in the magnetosphere, i.e. the exterior field is 
purely poloidal and is generated by the currents in the interior (type 
III). This assumption is justified by the estimate of Beloborodov 
(2009) that any twist in the magnetosphere should be dissipated on 
a timescale of 1 year or shorter. Under this constraint a toroidal 
magnetic field can only exist in regions where magnetic field lines 
close inside the star (type I), and vanishes otherwise (type II). We 
note that some of the work on magnetars (Colaiuda et al. 2008; 
Lander et al. 2012; Colaiuda & Kokkotas 2011b) use magnetic field 
configurations with a non-zero surface current to mimic unknown 
magneto spheric currents. We discuss below the relevance of these 
currents and their possible influence in our simulations, (iii) There 
are no current sheets or strong gradients in the magnetar interior be- 
cause there is enough evolution time since the magnetar' s birth to 
dissipate these. Therefore, we favour smooth functions with weak 
gradients across the star, i.e. we keep only low-order terms in the 
free functions of our magnetic field configurations. 

According to these restrictions we consider three magnetic 
field regions in magnetars: (i) the region in the star, where magnetic 



field lines close in the interior, can be of type I, II or III, (ii) the re- 
gion in the interior where magnetic field lines closing outside the 
star, can be of type II only (assuming that the currents generating 
the field originate in the interior of the star), and (iii) the exterior 
region of type III. It is convenient to treat all regions within the 
same framework using the current given by Eq. (21) for type I re- 
gions, and using the appropriate limits for and J'^ to construct 
the type II and III regions. 

We computed the equilibrium models numerically using the 
LORENE library^. To do so we have modified the MAGSTAR tool 
of LORENE to account for the general form of the current (21) in 
Ampere's law (16). This new routine is called MAGNETSTAR. We 
neglect the magnetic field terms in the space-time and hydrostatic 
equilibrium equations, and use the same equations and numerical 
methods as in Bocquet et al. (1995). We first solve the equation for 
the unmagnetized equilibrium configuration, and then solve Am- 
pere's law using the spherical spacetime and the density profile of 
the unmagnetized star as a background. We apply a fixed-point it- 
eration to solve Ampere's law, starting with an initial current ^i^i- 
The resulting magnetic field configuration depends on the functions 
M{A^p) and C^{A^), and on the initial current J*ni the fixed- 
point iteration. 

For M{A^p) we use the form 



M{A^) = 



l^ao A^ 
A 



-\-aiAi-\- a2 Al, 



+ac log 1 + 



+ de(A^ - A^{x)) 



(27) 



where A^ = A^jA^^^ , A^^^ is the maximum of |Ac^| inside the 
star, A^p (x) the value of A^ at the equator at a distance x from the 
center, 9 the Heaviside step function, and ao, ai, a2, ac, and d are 
constant coefficients, which control the amplitude of the perturba- 
tion M, and c a parameter to control the shape of the current. The 
field strength of the respective configuration is controlled by Mmax- 
Note that the expression following d leads to a delta function in the 
current, and corresponds to a infinitively thin circular current loop 
(Jackson 1998). For this case there exists an analytic solution which 
we use to avoid dealing with delta functions numerically. 
For Hcp{Acp) we use the form 



HU^^) = boA^ + biAl + b2Al , 



(28) 



where A^ 



■Al 



)/A™^"- with A-^^rf.max. 



maximum value of A^ at the surface. The resulting is zero at 
the surface, and thus it is confined inside the star. The case n = 1 
corresponds to one of the configurations in Lander & Jones (2012). 

Regarding the initial current J^- we consider three configura- 
tions 



spherical: J."^. = phjo , 
dipole-like: ^^^^ = phjor sin , 
quadrupole-like: J^.^. = phjor^ sin cos , 

where jo is a constant parameter of the same order of magnitude 
as the coefficients used for the current (ao, ai ...). The resulting 
magnetic field does not depend on the value of jo itself, and its 
strength is determined by the coefficients of Eqs. (27) and (28), i.e. 
ai,hi,ac, and d, respectively. 

Both the spherical and the dipole-like initial currents have 
even symmetry with respect to the equator. This also holds when 
they are rescaled to the physically relevant J{ip) — \fj\^J^ - This 



http://www.lorene.obspm.fr 



© 0000 RAS, MNRAS 000, 000-000 



Magneto-elastic oscillations of neutron stars II: exploring different magnetic field configurations 



J((^) has a maximum at the equator and resembles a current loop. 
The corresponding magnetic field configurations are similar to that 
of a simple current loop (compare with the models in Sec. 3.1.4). 
Since the current loop is an idealization giving rise to a pure dipole 
field in the far field, we characterize the magnetic field configura- 
tions generated by the spherical or dipole-like currents as dipole- 
like. We note that in the near field regime, where the currents have 
finite spatial extension, other multipoles are always present. Both 
the spherical and the dipole-like initial current lead to the same final 
equilibrium configurations. Similarly, we call the magnetic field of 
a quadrupole-like initial current, which has odd symmetry with re- 
spect to the equator, a quadrupole-like magnetic field. In this case 
the dipole moment vanishes and the far field is dominated by the 
quadrupole component. To obtain these quadrupole-like configura- 
tions one must use those currents that are an odd function of A^. 
Otherwise the fixed-point iteration converges to the same dipole- 
like magnetic field as in the previous case. We denote by D and Q 
the configurations where the converged magnetic field is predomi- 
nantly dipole-like or quadrupole-like, respectively. 

Note that for the magnetic field configurations with non- 
vanishing parameters ao, ai, and bo terms appear in the current 
that are linear or constant in A^. Therefore, solutions correspond- 
ing to different parameters ao, ai, and 60, can be linearly superim- 
posed in the current to obtain further solutions. Using this proce- 
dure it is possible to compute configurations with a mixed dipole- 
quadrupole-like magnetic field (see Section 4). 

We computed most of our initial models using the LORENE 
library and the procedure described above. However, in some cases 
where analytical solutions are known, we use these directly: the 
magnetic field originating from a ring current at a prescribed radius 
(see Section 3. 1.4), and force-free configurations confined to the 
crust according to Aguilera et al. (2008) (see Section 5). 



3 DIPOLE-LIKE CONFIGURATIONS 

In this section we study the effects of different dipole-like mag- 
netic field configurations on the Alfven spectrum of a magnetar. 
We consider either a purely poloidal component (Section 3.1) or a 
mixed poloidal-toroidal field whose poloidal component is dipolar- 
like (Section 3.2). Our analysis is mainly based on the application 
of the semi-analytic model to investigate the behaviour of purely 
toroidal oscillations. Using this model we compute the Alfven fre- 
quencies of individual field lines. The spectrum of the star is deter- 
mined by the set of the frequencies of all field lines and the corre- 
sponding overtones given as a function of position inside the star. 
The position can be described by the point x at which a field line 
crosses the equatorial plane or by the polar angle 6 along a given 
radius. We also employ some numerical simulations to evaluate the 
quality of our approach. 

3.1 Purely poloidal fields (type II) 

We first discuss magnetar configurations with a purely toroidal cur- 
rent (type II). We generate the initial models using the procedure 
described in Sections 2.2.1 and 2.2.2, and the LORENE library. We 
focus on type II magnetic fields, which are completely defined by 
Eq. (27). The corresponding model parameters are given in Table 1. 
We obtain the desired magnetic field strength by rescaling Mmax- 
In the first three models (Ao, Ai, and Cc) only one parameter in 
Eq. (27) is different from zero (ao, ai, or ac respectively). Addi- 
tionally, we consider two different combinations of these param- 
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Model 


ao 


ai 


a2 


ttc C 


Radius of 
ring in km 


Ao 


1 













Cc 











1 c 




Al 





1 
















0.43 


0.57 







F 


0.79 





0.21 








Rt^ - - - - - ^ 



Table 1. Values of the current parameters of the employed models according 
to Eq. (27). The models represent ring currents, denoting the radius 
of the ring in km. 



eters which are physically motivated. One of these, model F, re- 
sults in a very flat spectrum and is expected to produce long-lasting 
QPOs and large gaps between the spectra of adjacent overtones of 
the Alfven oscillations. Another model, model O, is chosen such 
that its spectrum reproduces some of the observed frequency ratios 
of SGR 1900+14. 

The currents J"^ and the magnetic field configurations of the 
different models are displayed in Fig. 1 in the top and bottom pan- 
els, respectively. The different currents lead to very similar mag- 
netic field configurations, all having their maximum either at the 
center of the star or along the equator, and a minimum in the center 
of the closed field lines. While the magnetic field lines look qualita- 
tively very similar in all cases, the current distribution can be quite 
different. 

Model Ao has a spherically symmetric current distribution. 
The corresponding magnetic field is stronger along the polar axis 
than along the equator. It can be considered as the limiting config- 
uration of model Cc for c ^ 00. Varying the constant c, the spatial 
distribution of the current changes from being spherically symmet- 
ric, c 1, to being ahgned with the polar axis, c <C 1 (see second 
and third panels in the top row of Fig. 1). The corresponding mag- 
netic fields behave similarly, i.e. they align with the polar axis for 
c <C 1 (bottom row). 

By choosing a non-spherical current distribution, model Ai 
(fourth column in Fig. 1), the maximum current is found at the 
equator near the outermost closed magnetic field line (x ^ 6 km). 
For this particular configuration the current vanishes at the polar 
axis. The resulting magnetic field is strongest in a central region 
oriented near the equator. Model O (fifth column in Fig. 1) looks 
very similar, but the current is somewhat more concentrated to- 
wards the equator, which leads to a weaker magnetic field close 
to the polar axis. Finally, model F has a flat spectrum (see below) 
and is characterized by two current maxima, one at the center and 
the other one at the equator at x 6 km. 

Due to their similar field configurations, the spectra of the cor- 
responding Alfven oscillations are quite similar (see Fig. 2). Near 
the polar axis (x = km) all spectra have a turning point (Ul) that 
can be either a maximum (Ao, Cc) or a minimum (Ai, F, O). Here, 
X is the radius where the magnetic field line crosses the equatorial 
plane. The presence of turning points in the spectra is of particular 
interest, since QPOs are expected at the corresponding frequencies 
(Levin 2007; Sotani et al. 2008; Cerda-Duran et al. 2009; Gabler 
et al. 2011, 2012). The QPOs close to the polar axis are called Up- 
per QPOs. 
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Figure 1. Upper panels: Current distribution J'-p for the model configurations (from left to right) Aq, Cio,Co.i, Ai, O, and F . Lower panels: Corresponding 
magnetic field configurations with selected field lines. The colour scale ranges from blue (minimum) to red (maximum) and represents the current strength (top 
panels) and the strength of the poloidal magnetic field (bottom panels), respectively. Here and henceforth we use cylindrical coordinates z and vj to indicate 
the height and radial distance from the z-axis, respectively. 



3.1.1 Single Upper turning point (Ul) 

Model Ao, black line in Fig. 2, was studied extensively in Gabler 
et al. (201 1, 2012) and has a second turning point at the closed field 
lines at X 6.5 km, which causes the appearance of the Lower 
QPOs (Cerda-Duran et al. 2009). The spectra of the two models 
Co.i and Cio, magenta and light blue lines in Fig. 2, are qualita- 
tively similar. The main difference between these two models is the 
gradient of the spectrum, which is stronger for the smaller value 
of c, because of the difference rates at which the magnetic field 
strength decreases with increasing x- For models where the elec- 
tric current is concentrated more along the polar axis, the magnetic 
field close to the axis is stronger than for spherically symmetric 
currents, and it decreases faster with increasing x- The stronger the 
magnetic field is close to the polar axis, the faster is the Alfven 
speed in this region, i.e. the frequencies near the polar axis in- 
crease, too. (cf. model Co.i). An abrupt change of the frequen- 
cies of the oscillations along neighbouring field lines is expected 
to lead to faster phase-mixing, which contributes to shorter lived 
QPOs (Levin 2007). Observing differences in the lifetimes of QPOs 
could in principle allow one to draw conclusions about the structure 
of the magnetic field close to the polar axis. Unfortunately, due to 
the numerical dissipation of our code we are not able to provide 
reliable estimates for the differences of the QPO lifetime caused by 
the different behaviour of the spectrum close to the polar axis. 

The spectra for models Cc have a minimum in the closed field 
line region (x ^ 5 km) only for c > 1 . Therefore, configurations 
with c < 1 do not have lower QPOs as observed in Gabler et al. 
(2012). However, this may not have strong observational conse- 
quences, because closed field lines couple only indirectly through 
the crust to the magnetosphere, i.e. we do not expect a strong mod- 
ulation of the emission occurring in the magnetosphere at the fre- 
quencies of the possible lower QPOs. 

The spectra of models Ai and O are qualitatively different (red 
and blue line in Fig. 2). The turning point at the polar axis is a min- 
imum for both models, and there exists a second turning point at 
X 5 km, i.e. we expect a new family of QPOs (U2) besides the 
previously observed Upper (U), Lower (L) and Edge (E) QPOs. 




2 4 6 

Crossing radius with equator % [km] 

Figure 2. Alfven spectra of the different magnetic field configurations 
given in Table 1 . The spectra are obtained with our semi-analytic model 
at a surface-averaged magnetic field strength of IQ-*^^ G. The spectra for 
X < 5 km correspond to open field lines, while those for x ^ 5 km to 
closed ones. Upper turning points (Ul and U2) are indicated black dots. The 
X-axis gives the point x where the field line crosses the equatorial plane. 



Edge QPOs are shorter-lived than turning-point QPOs. They are 
related to the edge of the spectrum at x ^ 5 km, where the outer- 
most open field line crosses the equator (see Gabler et al. 2012, for 
details). 

Sotani et al. (2008); Cerda-Duran et al. (2009) and Gabler 
et al. (2012) interpreted the observed QPO frequencies of 30, 92, 
and 150 Hz in SGR 1806-20 as the fundamental, the second and 
the fourth Alfven overtones of the Upper QPOs, respectively. Sim- 
ilar arguments hold for SGR 1900-1-14 and the observed frequen- 
cies of 28, 84, and 155Hz. The surface-averaged magnetic field 
strength that matches the fundamental frequencies of both SGRs 
for the chosen current configurations is listed in Table 2. It shows 
that for models where the current has a maximum close to the polar 
axis (Co.i, Cio, and Aq), the surface magnetic field strength neces- 
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Model 


turning point polar axis (Ul) 


second turning point (U2) 




-'-'lb 




Hi K 

-'-^15 


Hi 

-'-^15 




fn — 28 Hz 


fn — 30 Hz 


fn — 28 Hz 


fn — 30 Hz 


Ao 


4.9 


5.3 






Co.i 


3.4 


3.6 






Cio 


4.9 


5.3 






Ai 


6.8 


7.4 


4.7 


5.0 


O 


7.8 


8.3 


4.1 


4.3 


F 


5.5 


5.9 






R4 


4.1 


4.4 


3.4 


3.6 


R5 


5.4 


5.9 


3.0 


3.2 


Rq 


7.4 


7.9 






R7 


10.3 


11.0 








14.3 


15.2 







Table 2. Surface-averaged magnetic field strength needed to match a turn- 
ing point of the spectrum to /o = 28 Hz and /o = 30 Hz, respectively. 
The second and third column match the frequency of the turning point at 
the polar axis (Ul), while the last two columns match the frequency of the 
second turning point (U2) in the spectrum, if present (Ai, O, R4, and R^). 



80 



60 r 



40 



[1^ 



20- 



65.7 
52.7 
39.3 
24.8 



"0 2 4 6 8 

Crossing radius with equator % [km] 

Figure 3. Alfven spectrum of the magnetic field configuration F which leads 
to a very flat spectrum at B15 = 3.7. Red dashed lines indicate the frequen- 
cies of the crustal shear modes. The x-axis gives the point x where the field 
line crosses the equatorial plane. 



sary to match the observed 30 Hz in the SGR 1806-20 giant flare is 
lower than for the other models (F, Ai, and O). 



3.1.2 Two Upper turning points (Ul and U2) 

The presence of the second turning point in models Ai and O 
allows for an other interpretation of the observed QPOs of SGR 
1900+14: One could identify the fundamental QPO (Ulo) and its 
second overtone (UI2) with the observed frequencies of 28 Hz and 
84 Hz, whereas the observed frequencies of 53 Hz and 155 Hz could 
correspond to the fundamental QPO at the second turning point 
(U2o) and its second overtone (U22) (located at x 5 km)^. Model 
O was constructed such as to fulfil this relation between the ob- 
served frequencies. The corresponding model parameters (Table 1) 
are not very peculiar, i.e. given the current theoretical uncertainties 
model O is as valid as any other model we could have constructed. 



3.1.3 Nearly flat spectra: no shear modes inside continuum gaps 

We choose the parameters of model F to obtain a very flat spec- 
trum which is expected to produce very long-lasting QPOs in a 
large volume of the star. This happens because for almost all open 
field lines the QPOs have a very similar frequency. Additionally, 
the gap between successive overtones in the spectrum is maxi- 
mized, which may allow for crustal shear modes surviving in these 
gaps (Colaiuda & Kokkotas 2011b; van Hoven & Levin 2012) in- 
stead of being absorbed into the Alfven continuum of the core as 
in Gabler et al. (2012). Hence, we have computed the spectra of 
the fundamental Alfven oscillation and its first three overtones at 
= 3.7. These spectra are shown in Fig. 3 together with the fre- 
quencies of the crustal shear modes (horizontal red dashed lines). 
The / = 2 and / = 5 shear modes with 24.8 Hz and 65.7 Hz, 
respectively, clearly fall between the frequencies of Alfven over- 
tones, while the frequencies for the / = 3 and / = 4 modes lie 
just on top of the maximum or at the edge. To search for the crustal 




40 60 
Frequency [Hz] 



Figure 4. Fourier analysis of the osciUations at B15 = 3.7, 1.85, and 0.93. 
Magenta vertical lines indicate the frequencies of the first four crustal shear 
modes, whose precise values are given in Fig. 3. 



shear modes which may have been survived, we performed simula- 
tions covering a period of 0.25 s, calculated the overlap integrals^ 
with the pure crustal shear modes, and Fourier- analyzed them. The 
resulting Fourier transform for the / = 2 crustal shear mode at 
Bi5 = 3.7, 1.85, and 0.93 is given in Fig. 4. At B13 = 3.7 we find 
strong signals at the frequencies / ~ 20, 60, and 100 Hz. These 
are the frequencies expected from the semi-analytic model for the 
fundamental QPO and its second and fourth overtone (see Fig. 3). 
The frequencies of the overtones obey integer relations and do not 
behave as those of crustal shear modes which are proportional to 
fi ^ y^(/ — 1)(/ + 2)/o. None of the frequencies of the first three 
shear modes indicated by the magenta vertical lines in Fig. 4 coin- 
cides with an observed frequency of a QPO at B15 = 3.7. Fig. 4 
also shows that all observed QPOs scale directly with the magnetic 
field: the frequencies at B15 = 3.7 are twice as large as the cor- 
responding frequencies at B15 = 1.85 and are four times as large 



^ We use subscripts to indicate overtones of the fundamental QPO fre- 
quency. 



^ Overlap integrals are global measures of how strong a given shear mode 
is excited. For details we refer to Gabler et al. (2012). 
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G3 [km] 03 [km] 03 [km] 



Figure 5. Magnetic field configuration for a ring current at a given radius 
as indicated in the panels. The colour scale ranges from blue (minimum) to 
red (maximum), and gives the logarithm of the magnetic field strength. 




Crossing radius with equator % [km] 



Figure 6. Alfven spectra of magnetic field configurations obtained from 
ring currents R-oj at different radii vj. The spectra are calculated with 
our semi-analytic model for a surface-averaged magnetic field strength of 
Bi5 = 1. The X-axis gives the point x where the field line crosses the 
equatorial plane. 

as at 5i5 = 0.93. The matching of the frequency of the Alfven 
overtone at / 25 Hz and at 5i5 = 0.93 with that of the / = 2 
shear mode is pure coincidence. 

Qualitatively the same result holds when analyzing the over- 
lap integral for the / = 3 crustal shear mode. In this case all QPOs 
where the velocity field is symmetric with respect to the equato- 
rial plane are discovered and, again, no signal is present at the 
crustal shear modes. The analysis with the overlap integrals is lim- 
ited to oscillations having the same equatorial symmetry as the cor- 
responding shear mode, i.e. the overlap integral is zero for the / = 2 
(/ = 3) mode and an oscillation of odd (even) symmetry with re- 
spect to the equatorial plane. 

3.1.4 Ring current 

Next we consider the extreme case of a magnetic field generated by 
a circular current loop inside the star, whose analytical Newtonian 
solution can be found in Jackson (1998). A ring current can be seen 
as a limiting case of model Ai where the region of the current is 
shrunk to a ring. Placing the current at different radii we can con- 
struct different magnetic field configurations inside the neutron star 
(Fig. 5). 

The effects on the Alfven oscillations are illustrated in Fig. 6, 
which displays the spectra for models where the ring currents 
are located at t<7 = 4, 5, 6, 7, and 8 km from the center. In all cases 



Model 6o magnetic field strength [10^^ G] 

surface-averaged max. poloidal max. toroidal 



Bo 





3.5 


19 





Bi 


1 


3.6 


20 


3.5 


B2 


2 


4.0 


21 


7.1 


Bs 


5 


5.6 


24 


18 


Bio 


10 


7.6 


27 


33 


B20 


20 


8.9 


31 


54 



Table 3. Surface-averaged, maximum poloidal and toroidal magnetic field 
strength for mixed poloidal-toroidal configurations with parameters ao = 1 
and 60 as indicated. 




03 [km] 03 [km] 03 [km] 



Figure 7. Magnetic field lines of different mixed poloidal-toroidal models. 
The colour scale gives the strength of the toroidal magnetic field and ranges 
from white-blue (minimum) to red-black (maximum). The field of model 
Bo is purely poloidal. 

the spectra have a minimum at the polar axis x — ^ km (Ul), and 
the frequency of the open field lines increases with increasing x- 
Similar to model Ai there is a second maximum at x ^ 2.0km 
and X ^ 3.0 km for models i?4 and i?5, respectively. This sec- 
ond turning point may give rise to a new family of QPOs (U2) as 
discussed for models Ai and O in the previous subsection. 

The spectra associated with the closed field lines (right part 
of Fig. 6) show a strongly increasing frequency with increasing x 
until the location of the ring current is reached. There are additional 
turning points for models Ra, R5, and Rq ai x 2.5 km, x 
3.5 km and x 4.5km, respectively. However, the steepness of 
the spectra limits the existence of long-lived lower QPOs for these 
models, and for models R7, and Rg we do not expect lower QPOs 
at all. 



3.2 Mixed poloidal-toroidal field (type I) 

We turn now to models that have a mixed poloidal-toroidal field. 
In this case, the system of equations presented in Sec. 2.1 is no 
longer valid, as the toroidal component will couple axial and polar 
oscillations. However, since the toroidal component of our models 
is limited to the region of closed field lines inside the star, axial 
and polar oscillations still decouple at the linear level in the region 
of the open field lines. The coupling in the closed field line region 
will excite polar oscillations in the whole star. Fortunately, these do 
not influence the continuum at the linear level and we can still use 
our semi-analytic model to estimate the frequencies of the Alfven 
oscillation in the region of open field lines. 

Our models are obtained setting ao = 1 in Eq. (27) and 
Eq. (28) with different values for bo. The surface- averaged mag- 
netic field strengths and the maximum poloidal and toroidal field 
strength for the respective configurations are given in Table 3. The 
two components have similar strength at 60 ^ 10. The field lines 
and the toroidal magnetic field strength of selected models are plot- 
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pureQ Q/D=10 Q/D = 1 Q/D = 0.1 




02468 10 2468 10 2468 10 2468 10 



W [km] CD [km] 03 [km] 05 [km] 



2 4 6 8 

Crossing radius with equator % [km] 

Figure 8. Alfven spectra of torsional oscillations along the open field lines 
for mixed poloidal-toroidal models obtained with our semi-analytic 
method. See Table 3 for model details. The x-axis gives the point x where 
the field line crosses the equatorial plane. 



ted in Fig. 7, which shows that the area occupied by the closed field 
lines is largest in the case of zero toroidal field (model Bo). With 
increasing toroidal field this region shrinks and shifts towards the 
surface, see e.g. Lander et al. (2012). Additionally, the strength of 
the poloidal component increases at the closed field lines and at an 
limited region of the open field lines around the closed ones. 

In Fig. 8 we plot the spectra of the open field lines of our mixed 
field models. For zero toroidal field the Alfven spectrum is simi- 
lar to the black line in Fig. 2, but scaled at different frequencies. 
There is a maximum at the polar axis (x = km) and an edge at 
X 5 km. With increasing values of bo not only the toroidal field 
increases but also the poloidal one. This increase is stronger near 
the closed field lines, such that the frequencies at large x ('^ 5 km) 
increase faster than closer to the polar axis, leading to a flatter spec- 
trum in the case of model B2 (green line in the figure). With an 
even stronger toroidal component, the turning point at the polar 
axis changes from a maximum to a minimum, and the Alfven os- 
cillations near the closed field lines (x ^ 6 and 7 km for mod- 
els B5 and B20, respectively) show the highest frequencies. The 
larger bo is, the steeper the gradient of the spectrum becomes at the 
edge. In this case possible QPOs at the edge would last for shorter 
timescales than for model Bo where the gradient of the spectrum is 
moderate. 



4 MIXED DIPOLE-QUADRUPOLE-LIKE POLOIDAL 
CONFIGURATIONS 

To obtain configurations with a quadrupole-like contribution we 
only consider for simplicity purely poloidal magnetic fields, i.e. 
type II magnetic fields, by setting H{A^) equal to zero. Concern- 
ing the choice of parameters for the function M{A^), only config- 
urations that allow for a change of sign of J"^ across the equator 
and the corresponding initial data can be considered. In this section 
we choose ai / 0, while all other parameters are set to zero. 

As we have mentioned in Section 2.2.2 we construct mixed 
dipole-quadrupole-like magnetic field configurations by linearly 
superimposing both currents, which are obtained independently. 
For the purely quadrupole-like (dipole-like) configuration, Q (D), 
ai = 18.1 (ao = 11.66) and the averaged surface magnetic field 



Figure 9. Field lines of different mixed dipole-quadrupole-like (D-Q) con- 
figurations, the left panel showing the purely quadrupole-like case. The title 
of each panel indicates the Q-to-D component ratio. The red, dashed lines 
indicate the path along which the spectra of Fig. 10 have been obtained. 



strength is B = 6.0 x 10^^ G(B = 4.4 x 10^^ G). When construct- 
ing the mixed configurations we keep the dipole-like field strength 
fixed and add a quadrupole-like component rescaled by Q/D times 
the reference field. 

The field lines of selected mixed configurations are displayed 
in Fig. 9. The configuration with Q/D — 10 is predominantly 
quadrupole-like, while Q/D = 0.1 gives rise to a slightly de- 
formed dipole-like configuration. When the amplitudes of both 
components (Q/D = 1) are comparable, the field is very distorted, 
and has no symmetry with respect to the equatorial plane. 

We emphasize that we study the mixed dipole-quadrupole-like 
configurations without including a shear modulus in the crust, in 
order to isolate the effect that the inclusion of a quadrupole-like 
component has on the continuum of torsional oscillations. 



4.1 Alfven spectra using the semi-analytic model 

In Fig. 10 we plot the spectra of different combinations of Q and D. 
The X-axis gives the polar angle at which the field line crosses a 
circle with radius 7.5 km indicated with red, dashed lines in Fig. 9. 
This choice is motivated by capturing all field lines of the spec- 
trum in one single parametrization . Fortunately, the center of the 
region(s) of closed field lines lies almost on this circle for all con- 
figurations studied here. 

For Q/D — 0.1 (black lines in Fig. 10) the spectrum looks 
similar to the pure dipole-like case Ao . Qualitatively, we can iden- 
tify the end of the line at x 5.2 km in Fig. 2 with ^ ~ 1.0 rad 
in Fig. 10. Similarly, the local maximum at x 7.0 km in one fig- 
ure can be identified with the local maximum ^ 1.5 rad in the 
other figure. We note the asymmetry between the two hemispheres 
(0 < 7v/2 and > 7r/2) due to the presence of the quadrupole-like 
component. There is a maximum at the polar axis in the northern 
and the southern hemispheres for the open field lines (solid lines 
in the figure) and a minimum inside the region of closed field lines 
(dashed lines in the figure) at ^ ^ 1.1. The latter field lines also 
form the minimum at ^ ^ 1.7. Additionally, there is another max- 
imum at ^ ^ 1.5 that is associated with the closed field line which 
collapses into one point. Since this line has vanishing length, we do 
not expect to see a QPO at this location. 

When the quadrupole-like component is half as strong as the 
dipole-like one, model Q/D — 0.5 (green lines in Fig. 10), there 
appear two new extrema in the spectrum of the open field lines at 
^ 0.4 (U2) and ^ 0.8 (U3),respectively, which could pro- 
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Figure 10. Alfven spectra of different mixed configurations. Different line styles indicate regions of different types of field lines corresponding to Fig. 9. Solid 
lines correspond to open field lines, dashed lines to closed ones and the dotted line corresponds to open field lines only existing in the southern hemisphere for 
Q/D = 1.0. Selected turning points are indicated with black dots and labelled with U for Upper QPOs and L for Lower QPOs, respectively. The subscripts 
N and S stand for the hemisphere where the QPO appears, the northern < < 7t/2 (N) and the southern 7r/2 < < tt (S). The field strength of the 
dipole-like component is 4.4 x 10^^ G, and the quadrupole-like component varies from 0.6 x 10^"^ G at Q/D = 0.1 to 6.0 x 10^^ G at Q/D = 10.0, 
respectively. The pure quadrupole-like configuration is shown for B = 1.2 x 10^^ G. 



duce QPOs. Additionally, the spectrum becomes fairly asymmetric 
and the part due to the closed field lines (dashed lines in the fig- 
ure) lies almost completely in the northern hemisphere < 7t/2. 
At equal amplitudes of both magnetic fields, model Q/D = 1.0 
(red lines in Fig. 10), the part of the spectrum due to the open field 
lines simplifies: the turning points at ^ = 0.0 and ^ ^ 0.5 dis- 
appear, and only the one maximum at ^ ^ 0.75 (UIa^) remains. 
The closed field lines part of the spectrum (dashed lines) in the 
northern hemisphere is now limited to the interval 0.8 < < 1.5., 
and for 2.0 < < 2.5 there appears a new family of closed field 
lines disconnected from the previous ones at smaller 6. The new 
closed lines have lower fundamental Alfven frequencies, because 
the magnetic field in this area of the quadrupole-like component 
cancels partly with that of the dipole-like component, i.e. the re- 
sulting field is weaker than in the northern hemisphere. Addition- 
ally, there is a new family of open field lines (dotted red line in the 
figure) surrounding the closed lines at ^ > 7r/2, which we only 
observe for this particular combination Q/D = 1.0. They do not 
have any turning point, i.e. we expect at most very weak QPOs at 
the edges of the spectrum. 

For increasing values of Q/D the frequencies increase for 
field lines with > 0.5, because of the high magnetic field strength 
of the quadrupole-like component in this region. For the model with 
Q/D = 10.0 (magenta lines in Fig. 10) the symmetric spectrum of 
a pure quadrupole-like configuration is approached (see the blue 
lines in Fig. 10). For Q/D = 10.0 there is one turning point for 
the open field lines at^ ^ 0.7 (UIa^) that is formed by the same 
field lines leading to another turning point at ^ ^ 1 .4. They are also 
related to the two turning points (Uls) at ^ ^ 1.7 and ^ 2.4, re- 
spectively, caused by the mirrored field lines in the opposite hemi- 



sphere. The form of the spectrum associated with the closed field 
lines almost does not change in the northern hemisphere 6 < 7t/2 
compared to Q/D = 1.0, while the one in the southern hemisphere 
becomes symmetric to the one in the northern hemisphere. With 
increasing Q/D ratio, also the absolute values of the frequencies 
increase. 



4.2 Alfven QPOs using simulations 

For the study of the time evolution of the Alfven oscillations of the 
mixed dipole-quadrupole-like configurations we use the GRMHD 
code described in Sec. 2.1. The Fourier transforms of the strongest 
QPOs during the simulations for the configurations Q/D = 0.5, 
1.0, and 10.0 are shown in the top, middle, and bottom rows of 
Fig. 11, respectively. We find oscillations at the frequencies pre- 
dicted by the spectrum obtained with the semi- analytic model. The 
particular frequencies and names of the QPOs displayed in Fig. 1 1 
are indicated in the upper right of each panel. We selected those 
QPOs which are best suited for the visualization considering dif- 
ferent effects. For example, the particular initial perturbation ap- 
plied to the equilibrium model may have excited some QPOs more 
than others. Some QPOs may have similar frequencies, i.e. their 
signals may overlap in the corresponding Fourier transform, which 
has a finite resolution. As expected, the spatial structure of the 
QPOs shown in the figure follows the structure of the magnetic field 
lines. For the Q/D = 0.5 (top row) the maximum amplitude of the 
Fourier transform for the UIq QPO (first column) is located near 
the polar axis, for U3i (fourth column) it is at the open field lines 
nearest to the closed field lines, while for U22 and U24 (second 
and fifth columns) the maximum amplitudes are found in between 
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Figure 11. Fourier transforms of the Alfven QPOs for configurations QjD = 0.5 (top), Q/D = 1 (middle), and Q/D = 10 (bottom), respectively. The 
plotted QPOs are: f/lo, f/22, Lq, f/3i, and [724 (top), Ulsi, Lso, UInq, L^q, and Ul^i (middle) UInq, Lso, Ulsi, L^q, and f/ljvi (bottom). The 
subscript indicates the number of the overtone, with indicating the fundamental QPO. In all panels the colour scale ranges from white-blue (minimum) to 
red-black (maximum) and gives the logarithm of the Fourier amplitude of the QPOs at the frequency denoted in the upper right of each panel. 



these two regions. The oscillation of Lq (third column) is limited 
to the closed field lines. In the third column we can also see three 
overtones of the different Upper QPOs: UI2 (three nodes near the 
polar axis), U23 (four nodes), and U3i (two nodes along the open 
field lines close to the closed ones), which all have approximately 
the same frequency. 



The QPOs of the configuration with Q/D = 1.0 are plotted 
in the middle row of Fig. 11 and are labelled Ulsi, Lso, UIno, 
Lno, and UIni (from left to right, respectively). The additional 
subscripts N and S indicate the northern or southern hemisphere 
where different QPOs are localized. Ulso and Lso are the funda- 
mental oscillations of the open and closed field lines in the southern 
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hemisphere, while UIno, Lno, and C/livi are limited to the north- 
em one. 

For Q/D = 10 (bottom row) both hemispheres approach a 
similar magnetic field configuration, which results in an almost 
symmetric distribution of the oscillations. In most of the panels 
(2-4) of this row we see a strong Fourier amplitude in both hemi- 
spheres. This is expected because the frequencies of the corre- 
sponding QPOs in the northern and southern hemispheres are al- 
most equal. However, there are still some small differences in the 
frequencies, e.g. between Ulsi (f ~ 8.1Hz) and UIni (f ~ 
8.9 Hz), or between Lso (f 7.7Hz) and Lno (f 8.3 Hz), 
respectively. 



5 MAGNETIC FIELD CONFIGURATIONS CONFINED 
TO THE CRUST 

We now turn to the study of magnetic field configurations that are 
confined to the crustal region of the neutron star. In such cases the 
crustal shear modes cannot be absorbed by the core. Our choice 
is motivated by the possible presence of superconducting protons 
in the core of the neutron star (Page et al. 2011; Shternin et al. 
2011) which could expel the magnetic flux due to the Meissner- 
Ochsenfeld effect. From theoretical studies (Baym et al. 1969) one 
would expect the protons to form a superconductor of type II, the 
magnetic field being able to penetrate the superconducting region 
in form of flux tubes. However, due to the uncertainty of the EoS 
at supranuclear densities, the opposite situation cannot be ruled 
out completely. In the present context only consider study the ef- 
fects of type I superconductivity on the Alfven spectrum of tor- 
sional modes. The only other existing study of such a configuration 
(but with an additional toroidal component and assuming a discrete 
spectrum of modes) was presented by Sotani et al. (2008). 

We follow Aguilera et al. (2008) who give a description of dif- 
ferent axisymmetric magnetic field configurations. We neglect the 
effects of general relativity on the magnetic field and the influence 
of the magnetic field on the neutron star structure. This should be 
a valid approach, in particular, for fields confined to the crust con- 
taining only a few per cent of the total stellar mass. Moreover, we 
assume that the magnetic field has relaxed to some equilibrium be- 
fore the crust crystallized. After crystallization the magnetic field 
is frozen into the crust. For simplicity, we consider purely poloidal 
fields only. These models can be described by the type III mag- 
netic field of Section 2.2.1. The corresponding currents outside the 
crust determining the magnetic field structure are the superconduct- 
ing currents at the crust-core interface. Details on how to construct 
magnetic fields confined to the crust are given in Appendix A. 

We study different configurations, namely magnetic fields 
matched to an exterior dipole, quadrupole, or octupole field. The 
amplitudes of the magnetic field at the surface of such configura- 
tions a have different angular dependence. To compare the results, 
we thus label the different configurations with their averaged mag- 
netic field strength at the surface of the star, and use this value as 
reference magnetic field strength. 

The field topology of the configurations and the corresponding 
spectra of the purely Alfven oscillations BiB = 10^^ G are given in 
Fig. 12. The parts of the Alfven spectrum associated with a particu- 
lar family of field lines share the same colour in the graphs showing 
the spectrum and the magnetic field configuration. If the magnetic 
field is confined to the crust the maximum field strength is found 
near the crust-core interface. The maximum strength for the dipole- 
like (quadrupole-like, octupole-like) configuration is about 13 (8, 



8) times the average value at the surface. This is a factor of a few 
stronger than the crustal field of configurations penetrating the core. 
Consequently, the structure of the magnetic field varies on smaller 
scales. There are field lines which close inside the crust, e.g. near 
the equator for the dipole-like (red Unes in Fig. 12) and octupole- 
like (green lines) configurations, or at ^ ~ 60° for quadrupole-like 
and octupole-like configurations (see blue lines in Fig. 12). Because 
of stronger magnetic fields and stronger field gradients on smaller 
spatial scales, we need higher numerical grid resolution for our sim- 
ulations (see below) and expect higher frequencies of the Alfven 
oscillations. 

In the following we will refer to magnetic field configurations 
matched to an exterior dipole as configuration D, to an exterior 
quadrupole as Q, and to an octupole as O, respectively. 



5.1 Alfven spectra with the semi-analytic model (without 
shear) 

Let us first discuss the spectra obtained with the semi-analytic 
model of Cerda-Duran et al. (2009) and Gabler et al. (2012) dis- 
played in the lower panels of Fig. 12. First, we note that for con- 
figurations matched to an exterior quadrupole and octupole, there 
are no turning points present inside the continuum of the open field 
lines. Therefore, one would expect QPOs only at the edges of the 
different branches of the continua. For the dipole-like case we ob- 
serve a turning point only at the very high frequencies at field lines 
located around the closed field lines at the equator (corresponding 
to the upper end of the red coloured part of the spectrum). Second, 
the frequencies of the fundamental oscillation of the short closed 
field fines (> 100 Hz) are higher than the frequencies we observe 
for a magnetic field penetrating the core (a few Hz at 10^^ G). The 
frequencies of the field lines near the polar axis are similar to those 
obtained for a magnetic field configuration penetrating the core. 
Note that there is a edge at ^ = 0, and that the corresponding fre- 
quencies are finite. 



5.2 Alfven QPOs from simulations (without shear) 

The semi-analytic model shows that the spectra of magnetic fields 
confined to the crust are very steep near the edges and posses a 
turning point for configuration D at very high frequencies. More- 
over, the field lines near the edges are packed very densely along 
the crust-core interface leading to enhanced phase mixing. There- 
fore, we expect more difficulties to observe QPOs with significant 
amplitudes in our numerical simulations. As mentioned before, the 
requirements on the grid resolution are more stringent, because of 
the finer spatial structure of the magnetic field. The grid resolution 
employed for the time evolution of the three magnetic field config- 
urations (at B = 10^^ G) is 30 x 80 (r x 0) for [rcc, n] x [0, tt]. 
All models have been evolved for up to 5 s. 

Indeed, the Fourier analysis displayed in Fig. 13 shows less 
significant amphtudes of the QPOs than in the case of global mag- 
netic field configurations penetrating the core. Nevertheless, we can 
identify some QPOs as local maxima of the Fourier amplitude. The 
QPOs for the dipole-like configuration (upper panels) have their 
maximum amplitudes near the polar axis. This is expected, because 
the edge of the continuum is located there. With the numerical reso- 
lution used we were not able to see the other edge of the continuum 
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Figure 12. Upper panels: Magnetic field configurations for models matched to an exterior dipole (left), quadrupole (middle), and octupole (right), respectively. 
Lower panels: Alfven spectra of the magnetic field configurations given in the corresponding upper panel obtained with the semi-analytic model for an averaged 
surface magnetic field ofB ~ 10^^ G. Corresponding field lines in the upper and lower panels have the same colour. 



at much higher frequencies, because the associated field lines near 
the equator have a very small spatial extension^ . 

The quadrupole-like magnetic field has closed field lines near 
^ ~ 1.0 rad (see Fig. 12). The corresponding frequencies of the 
continuum are > 100 Hz. Unfortunately, the affordable grid res- 
olution is too low to analyze potential QPOs for these field lines. 
Therefore, we find QPOs only for the continuum of the open field 
lines. These are best seen in the plot for 4.3 Hz in the middle panel 
of Fig. 13. Only open field lines oscillate at the given frequency 
and have large Fourier amplitudes near the polar axis and equator. 
Closed field lines, on the other hand, do not participate in these os- 
cillations, as the Fourier amplitude vanishes in the corresponding 
regions. 

When the magnetic field is matched to an exterior octupole, 
we find oscillations of the open field lines (best visible for / = 
3.9 Hz in the bottom panel of Fig. 13). The red-coloured field lines 
in the right panel of Fig. 12 do not extend to the equator, and the 
corresponding QPO has a large amplitude around > 1.0 rad. 
At 0.9 Hz we find oscillations with significant amplitudes also for 
the field lines near the equator. At almost all frequencies there is 
a non-vanishing Fourier amplitude at the closed field lines around 
^ 0.75 and = 7r/2, which is probably caused by the too low 
grid resolution. 

In none of the cases considered is it possible to clearly asso- 
ciate the observed QPOs with features of the semi-analytic spec- 
trum obtained in Sec. 5.1, because the frequencies do not match 



^ To resolve these smaller scales, the grid resolution has to be refined sig- 
nificantly. However, we need to evolve the models for several dynamical 
time scales of the low frequency oscillations around 1 Hz, which takes the 
order of several days of computing time with the current resolution. Since 
we are mainly interested in the frequency range around 30 Hz, we had to 
find a compromise between the computational effort and the scanned fre- 
quency range. 



perfectly. Nevertheless, a conservative interpretation would be that 
some QPOs which are associated with the open field lines near the 
polar axis correspond to the edges of the continuum and their over- 
tones. The closed field lines are expected to oscillate at frequencies 
above 100 Hz which we cannot confirm because of lack of spatial 
resolution. Furthermore, we are mainly interested in the frequency 
range below 100 Hz and models where the elastic properties of 
the crust are included (see Section 5.3). Additionally, as a conse- 
quence of the finer spatial structure of the magnetic field confined 
to the crust, there is a much stronger numerical coupling of dif- 
ferent field lines, which leads to enhanced damping and weaker 
QPOs. We expect this weakness of our approach to be less impor- 
tant in the coupled scenario, because for sufficiently low magnetic 
fields, the oscillations will be dominated by the shear in the crust. 
The corresponding modes are sufficiently resolved at the numerical 
resolution applied here. 

5.3 Coupled magneto-elastic oscillations 

We now investigate the behaviour of coupled magneto-elastic os- 
cillations of the crust for the magnetic field configurations consid- 
ered above. Therefore, we performed a number of simulations with 
initial data consisting of a general perturbation with an angular de- 
pendence of the sum of the / = 2 and 3 vector spherical harmonics 
at different magnetic field strengths. In the following we will use 
the letter I' to label shear-mode-like QPOs inside the crust, where 
I' is the number of nodes in ^-direction. For non- vanishing mag- 
netic fields, these QPOs are no longer normal mode oscillations, 
but QPOs with time-varying spatial structure. Therefore, they do 
not have the angular dependence of individual vector spherical har- 
monics and they cannot be described by simple analytic functions. 
The chosen grid resolution for the simulations depends on the se- 
lected model and on the radius at which the crust-core interface 
Tec is located. The smallest resolution used consisted of 60 x 80 



© 0000 RAS, MNRAS 000, 000-000 



14 Michael Gabler, Pablo Cerdd-Durdn, Jose A. Font, Ewald Miiller and Nikolaos Stergioulas 



dipole-like D 




dipole quadrupole octupole 



02468 02468 02468 10 
05 [km] 05 [km] 05 [km] 



quadrupole-like Q 




02468 02468 02468 10 
05 [km] 05 [km] 05 [km] 

octupole-like O 




02468 02468 02468 10 
05 [km] 05 [km] 05 [km] 

Figure 13. Spatial structure of the low-frequency QPOs obtained by a 
Fourier analysis of different magnetic field configurations. Upper panels: 
configuration D; Middle panels: configuration Q: Lower panels: configu- 
ration O. The corresponding frequencies are given, too. The colour scale 
ranges from white (minimum) via red to black (maximum). 
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Figure 14. Fourier amplitude of the = 2 crustal shear mode and its 
magneto-elastic generalization for different magnetic field configurations 
and different field strengths. The colour scale ranges from white (minimum) 
via red to black (maximum).Bdi = 1.7 x lO^^ G, Bq^ad = 1.9x IQi'^ G, 
and Boctu = 1.2 x lO^^G, respectively. 



(r X 0) points for [rcc, ^s] x [0, tt]. The integration time is t 2 s 
2ii B — 10^^ G. For stronger fields, the integration time scales in- 
versely with the magnetic field strength. 

The results of the Fourier analysis of these simulations are 
plotted in Fig. 14, where we focus on the Fourier amplitude of 
the magneto-elastic generalization of the I' = 2 shear mode for 
different magnetic field strengths and different field configura- 
tions (dipole-like, quadrupole-like, and octupole-like). The aver- 
aged surface magnetic field strengths are ^di = 1.7 x 10^^ G, 
Bquad = 1.9 X 10^^ G, and Boctu = 1.2 x 10^^ G, respectively. 
When making general statements holding for all three magnetic 
field configurations we use the label B^. 



AtB — B^ (first row of Fig. 14) the structure of the magneto- 
elastic QPOs is very similar to that of the purely crustal shear 
modes for all three field configurations, and for all shear modes. 
The corresponding panels in Fig. 14 show how the purely crustal 
shear modes approximately look like. For larger magnetic field 
strengths, the assumed magnetic field configuration changes the 
structure of the shear modes. For example, in the presence of a 
dipole-like surface field, the I' = 2 QPO is compressed towards 
those field lines entering the star closer to the pole. The stronger 
the magnetic field is the closer to the polar axis the QPO are lo- 
cated (see at the leftmost panels of the figure from top to bottom). 
This resembles the behaviour of the QPOs of purely Alfven os- 
cillations for this configuration (see Fig. 13). On the other hand. 
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Figure 15. Dependence of the frequencies of the n = 0, 1' = {2, 3, 4, 5} 
QPOs on the magnetic field strength for different magnetic field config- 
urations which are confined to the crust. The purely crustal shear mode 
frequencies are given at the y-axis label as 24.6, 38.9, 52.2, and 65.1 Hz, 
respectively. The lines are fits according to Eq. (30), and the corresponding 
parameters are given in Table 4. 



the QPO structure of quadrupole-like surface fields (central pan- 
els) behaves in the opposite way, i.e. the maximum amplitude is 
shifted towards the equator with increasing magnetic field strength. 
Finally, for an octupole-like configuration (rightmost panels), the 
maximum Fourier amplitude remains at approximately its origi- 
nal position but it becomes narrower with increasing magnetic field 
strength. For the corresponding change in the frequency of crustal 
shear modes with I' > 2 (not shown in the figure) we observe a 
similar behaviour depending on the particular field configuration: 
The stronger the magnetic field, the more the structure of the QPO 
tends to resemble the structure of that field, and the shear modes 
are expelled from regions of closed field lines (see e.g. the region 
close to the equator for configurations D and O in Fig. 13). 

A common feature of all configurations is that with increas- 
ing magnetic field strength the identification of QPOs becomes 
more difficult because their amplitudes decrease. That is why some 
QPOs are not present at all field strengths. In particular, we have 
not succeeded in identifying the following QPOs: I' — 4 at 5 = 
10 X Bquad and B = 20 x 5quad, and = 5 at 5 = 20 x 5octu, 
respectively. For even stronger magnetic fields 5 = 50 x Bx it is 
almost impossible to identify the magneto-elastic generalizations of 
the crustal shear modes, as the evolution is completely dominated 
by the magnetic field. 

In order to understand how the frequencies of the shear modes 
change, when they become QPOs in the presence of the magnetic 
field we display in Fig. 15 the QPO frequencies as a function of 
the averaged surface magnetic field strength. Different colours in- 
dicate different QPOs associated with the zero-magnetic-field shear 
modes (black: I' = 2; red: I' = 3; magenta: I' = 4; blue: I' = 5). 
The line style indicates the magnetic field configuration (solid: 
dipole-like; dashed: quadrupole-like; dotted: octupole-like). At a 
surface magnetic field strength of about a few x 10^^ G the frequen- 
cies of the QPOs begin to deviate significantly from those of purely 
shear modes. For the = {2, 3, 4, 5} modes the zero-magnetic- 
field frequencies are 24.6, 38.9, 52.2, and 65.1Hz, respectively, 
and are marked along the 2/- axis 

Sotani et al. (2008) studied similar magnetic field configura- 
tions and found an increase of the frequency of the QPOs with in- 
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Table 4. Parameters Bcrit and 6^/ of Eq. (30) for different magnetic field 
configurations and different number of nodes in ^-direction V . 



creasing magnetic field. Sotani et al. (2007) give an approximate 
formula of the dependence of the frequency on the magnetic field 



fv{B) 



= \ 1 + ai, 



(29) 



where /° is the frequency of the purely shear eigenmode, and 
B/j^ = 4 X lO^^G. When fitting the frequencies obtained in our 
simulations with this formula we find a fairly inaccurate correspon- 
dence. Therefore, we generalized Eq. (29) to 



fAB) 



1 + 



B 



(30) 



where Bcrit, is the magnetic field strength at which the influence 
of the magnetic field becomes important. The parameters B^vit and 
bi' are given in Table 4, and the corresponding fits are displayed in 
Fig. 15. Note that we only used the values for B ^ 10 x Bx for the 
fits because either we could not find shear-mode-like oscillations 
for stronger fields or their identification was difficult. 

The observed exponents deviate significantly from hi' — 2, 
the majority being 4/3 < bii < 5/3. For very strong magnetic 
fields (about several x 10^^ G, see Gabler et al. 2012) one would 
expect a transition to Alfven oscillations, where the frequency de- 
pends linearly on i.e. 6^/ ^ 2. However, in the presence of such 
strong fields, shear oscillations no longer exist. In the regime stud- 
ied here we expect a significant influence of the shear modulus such 
that the asymptotic regime of 6^/ = 2 is not reached yet, which may 
explain the deviation from bi^ = 2. 

While there is no clear dependence of the exponent bi' on the 
magnetic field configurations, one should note that the critical sur- 
face field Bcrit is slightly lower for the dipole-like configuration 
than that for configurations Q and O for all QPOs and for all /^ 
This can be explained with the magnetic field structure. First, the 
configuration D has the strongest magnetic field inside the crust 
(^ 13 X Bsuvf) for a given surface field compared to the other two 
(^ 8 X Bsuvi)- Second, the n = shear modes^ travel predomi- 
nantly in ^-direction. Therefore, the ^-component of the magnetic 
field should be the dominant one to enhance the mode propagation 
velocity and, hence, to increase the frequency of the QPOs. In con- 
trast to configuration Q and O, which both have more than one node 



n is the number of nodes in radial direction. 
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of the ^-component of the magnetic field in ^-direction, the dipole- 
Uke configuration has only one node at the polar axis. Thus, the 
averaged absolute value of the ^-component is stronger for config- 
uration D and we expect it to influence the magneto-elastic waves 
more strongly than in case of the other two configurations. 

We can extrapolate our results given by Eq. (30) to different 
values of the shear modulus and test whether a change of this pa- 
rameter can explain the lowest observed frequencies of 18, 26, and 
30 Hz in SGR 1806-20, as proposed by Steiner & Watts (2009). 
If the shear modulus /Xs is divided by a factor F, the system of 
equations (5), (10), and (11) conserves its form if we rescale the 
background magnetic field hy the square root of this factor and 
simultaneously multiply the time by a/F: 



dt 

dt 
at 

where i t^F, t^, 



dt 



dx^ 



(31) 
(32) 
(33) 



B^/VF, and As 



IJis/F, respectively. 

The rescaling of the time coordinate implies an inverse rescal- 
ing of the frequency. Therefore, we obtain the expected result that 
in the absence of a magnetic field the frequencies of the crustal 
shear modes scale as 1/ a/F. A decrease of the shear modulus by a 
factor F leads to a corresponding decrease of the QPO frequency 
by a/F. Consequently, the critical magnetic field for an equilibrium 
model with a rescaled shear modulus (divided by F) has to be di- 
vided by \/F. Knowing 5crit and hi' for an equilibrium model we 
can compute the frequency of the QPOs for different shear moduli 
by 



fv{B) 



Vf 



1 + 



B 



Bcrit/VF 



(34) 



The corresponding curves for the dipole-like background configu- 
ration, the lowest QPO with I' = 2, and F = {2, 4, 8, 16, 100} are 
given in Fig. 16. Clearly, the decrease of the shear modulus leads 
to a decrease of Bcvn, i.e. the magnetic field becomes dynamically 
dominate at lower magnetic field strength than for large /is- With 
the dipole-like configuration we are not able to find oscillations at 
frequencies of 18 , 26 or 30 Hz at any value of the shear modulus 
up to l/16th of its tabulated value for B > 6 x 10^^ G, which is 
the lower bound for the dipole-like magnetic field of the magne- 
tars showing giant flares from spin-down estimates. Only for the 
unrealistically small value of jls = /Xs/lOO we reach the observed 
frequencies / = 30 and / = 28 Hz for a magnetic field strength of 
Bi3 < 1. 

To summarize the findings of this section, we see a signifi- 
cant increase in the frequency of the magneto-elastic generalization 
of the crustal shear modes with increasing magnetic field strength. 
The details of this increase, like the dependence on the magnetic 
field strength or the spatial structure of the resulting QPOs, de- 
pend on the particular magnetic field configuration. For higher mul- 
tipoles we find a weaker increase of the QPO frequencies with 
the magnetic field strength than for dipole-like fields (see Fig. 15). 
However, the estimate of the magnetic field strength with the spin- 
down formula gives a lower estimate on the dipole-like component 
of the magnetic field. Therefore, the results of this section can be 
seen as providing a lower estimate of the frequency increase if there 




Averaged surface magnetic field strength [10 G] 



Figure 16. Extrapolation of the QPO frequencies using Eq. (34) for dif- 
ferent values of the shear modulus /is in the presence of a dipole- 
like magnetic field for the QPO with V = 2. The red lines in- 
dicate the lowest observed frequencies of the QPOs for the range 
of magnetic field strengths obtained using spin-down estimates for 
the three magnetars showing giant flares (B15 = 0.6... 2.1, see 
http://www.physics.mcgill.ca/~pulsar/magnetar/main.html). 



are additional multipoles present in the magnetic field structure. At 
surface field strengths of B15 ^ 1, expected from the spin-down 
estimates, the frequencies of the crustal shear modes are shifted to 
such high values that it is difficult to explain the lowest frequencies 
of the QPOs observed in SGRs even for strongly reduced shear 
moduli. 



5.4 Different equilibrium models 

In addition to the APR-i-DH EoS used for all the results discussed 
so far in this paper, we have also investigated neutron star models 
of 1.4 solar masses obtained with the APR-hNV, L-i-DH and L-hNV 
EoS, respectively. For all three EoS we have performed simulations 
at magnetic field strengths of B15 = 0.1, 0.5, 1.0, 2.0, and 5.0, re- 
spectively. As discussed in the previous section, it was highly prob- 
lematic to identify the QPOs corresponding to the pure shear modes 
in the zero field limit for the strongest magnetic fields considered, 
because strong magnetic fields significantly change the QPOs. In 
particular those QPOs having a maximum at the equator (I' = 3 
and = 5) are difficult to detect for the dipole-like field config- 
uration. The Alfven oscillations along the closed field lines in this 
region have very different structure and frequencies, and, thus, sup- 
press shear-like oscillations. 

The corresponding parameters Bcrit and bi^ for the fitting 
functions defined in Eq. (30) for B15 < 1.5 are given in Table 5 for 
the set of EoS and magnetic field amplitudes. There is a clear, al- 
beit weak, dependence of the two parameters on the QPO structure 
characterized by /^ For low I' the critical magnetic field is weaker 
than for high l\ and exactly the opposite holds for the exponent 6^/ , 
which decreases with increasing l\ The influence of Bcrit on the 
crustal EoS is of the same order of magnitude as the difference be- 
tween different I'. The DH EoS leads to approximately 30% lower 
values of 5crit than the NV EoS. This strong effect of the mag- 
netic field for models with the DH EoS can be easily understood 
because the NV EoS leads to models with larger shear moduli. The 
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Table 5. Parameters B^vit and 6^/ of Eq. (30) for a magnetic field configu- 
ration matched to an exterior dipole for different EoS and different number 
of nodes in ^-direction V , respectively. 



influence of the core EoS is less pronounced and one recovers very 
similar 5crit and by for the two examples studied here. 

The dependence on the particular choice of the EoS is not very 
strong, and in all cases we find that the frequencies increase signif- 
icantly starting at a few 10^^ G. Therefore, for these configurations 
and all EoS analyzed, we cannot interpret the lowest observed QPO 
frequencies in terms of shear-like QPOs. Additionally, without the 
identification of the fundamental shear mode, it is not possible to 
explain the order of the higher frequency QPOs (/ > 30 Hz) with 
shear modes. If the = 2 mode has a higher frequency of the or- 
der of the frequencies of the higher V modes does not match with 
observations. Another possibility would be that the fundamental is 
not excited, but in this case one would have to explain why this 
happens. 



6 DISCUSSION 

We studied the effect of different magnetic field configurations on 
the magneto-elastic oscillations of magnetars. For this purpose we 
have constructed magnetized equilibrium configurations generaliz- 
ing the MAGSTAR routine of the LORENE library to include various 
descriptions of the current generating the magnetic field: MAGNET- 
STAR^. The oscillation spectrum, dominated by the continuum of 
the core, was studied by means of a semi-analytic model (Cerda- 
Duran et al. 2009; Gabler et al. 2012) and numerical GRMHD sim- 
ulations. 

Table 6 summarizes the main results for the magnetic field 
configurations considered in this paper. For each model we identify 
the fundamental Alfven QPO as the turning point in the Alfven con- 
tinuum at the lowest frequency. We identify this mode with some 
of the observed QPOs, the 30 Hz oscillation of SGR 1806-20 and 
28 Hz of SGR 1900+1, and we estimate the equivalent dipole mag- 
netic field strength of each configuration to match both frequencies. 
We define the equivalent dipole magnetic field, B, as the magnetic 

^ Our generalization of MAGSTAR is publicly available in the standard 
LORENE library. 



Model QPO Bf^^^ 



AO 


Ul 


2.1 


2.3 




Ul 


1.3 


1.4 


Cm 


Ul 


2.1 


2.2 


Al 


Ul 


3.1 


3.3 


Al 


U2 


2.1 


2.3 


O 


Ul 


3.6 


3.8 


O 


U2 


1.9 


2.0 


F 


Ul 


2.4 


2.5 


Ra 


Ul 


1.1 


1.2 


R4 


U2 


0.85 


0.93 


Rt^ 


Ul 


1.5 


1.6 


R5 


U2 


0.80 


0.85 


Re 


Ul 


2.0 


2.1 


R7 


Ul 


2.7 


2.8 


Rs 


Ul 


3.5 


3.8 


60 = 


Ul 


2.1 


2.3 


60 = 1 


Ul 


2.2 


2.4 


bo = 2 


Ul 


2.4 


2.6 


60 = 5 


Ul 


2.9 


3.1 


bo = 10 


Ul 


3.5 


3.8 


bo = 20 


Ul 


3.7 


4.0 


Q/D=0.1 


Ul 


2.1 


2.2 


Q/D=0.5 


Ul 


2.5 


2.7 


Q/D=0.5 


U2 


3.5 


3.7 


Q/D=0.5 


U3 


1.5 


1.6 


Q/D=1.0 


Uliv 


0.93 


0.99 


Q/D=1.0 


Vis 


4.0 


4.3 


Q/D=2.0 


Uliv 


0.53 


0.57 


Q/D=2.0 


Vis 


0.87 


0.93 


Q/D=5.0 


Uliv 


0.24 


0.25 


Q/D=5.0 


Vis 


0.28 


0.30 


Q/D=10.0 


Uliv 


0.12 


0.13 


Q/D=10.0 


UI5 


0.13 


0.14 


Q 


Ul 


0.076 


0.081 



Table 6. Equivalent field strength B to match the frequencies at / = 28 and 
/ = 30 Hz with the QPO indicated in the table. B is defined by m = BR^ 
for a uniformly magnetized sphere, where m is the magnetic moment and 
R is the radius of the star. The QPOs Unx are the Upper QPOs of the 
open field lines, n indicates the number of the corresponding turning points 
labelled from n = 1 at the polar axis to maximum n = 3 at larger 0, 
and X gives the division into northern (X = N) and southern (X = S) 
hemisphere, if there exist different QPOs in both hemispheres. 



field strength at the surface of a Newtonian, uniformly magnetized 
sphere having the same dipole magnetic moment^ m as the con- 
figuration we want to describe, which can be directly compared 
with the magnetic field estimates from spin-down measurements, 

through B[G] > 3.2 x 10^^ (^[s]P)^^^ (Lorimer & Kramer 
2004). Thus, the magnetic field estimates in Table 6 are directly 
comparable to estimates of the magnetic field from the spin down 
of the neutron star. 

We summarize and discuss our main results next: 

• For poloidal, purely dipole-like configurations generated by a 
toroidal current inside the star (type II) we observe that the Alfven 

^ The value of m is provided in the output of the LORENE library, and the 
corresponding magnetic field strength can be computed from m = BR^. 
For the ring current m is given by m = I A, where / is the current and A 
the area of the current loop. 
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continuum always contains turning points at the polar axis (Ul). 
Some models (Al, O, R4, and R5) have a second turning point 
(U2) in the region of open field lines inside the star. It seems very 
promising to identify some of the frequencies observed in the SGR 
giant flares, which are in the integer relation of 1 :3:5 to each other^, 
with the fundamental Alfven Upper QPO at the polar axis and its 
second and fourth overtone. As can be seen in Table 6, the required 
equivalent field strength is in the range of 0.8 to 3. 8x 10^^ G. This 
is consistent with the spin-down estimates of B15 > 2.4 for SGR 
1806-20 and Bi^ > 0.7 for SGR 1900+14 (where B15 is B in units 
of 10^^ G), respectively. 

• The appearance of more than one turning point in the spec- 
trum of some of the configurations allows us to relate additional 
observed frequencies to the torsional Alfven oscillations. Model O 
was constructed such that we can match all frequencies observed 
in SGR 1900-hl4. It is possible to identify the 28 Hz with Ul at the 
polar axis, and the 84 Hz= 3 x 28 Hz with its second overtone. 
Additionally, the 53 Hz QPO could be associated to U2 away from 
the polar axis at x 5 km (see blue line in Fig. 2) and 155 Hz 
~ 3 X 53 Hz could be its second overtone. Interpreting the 30 Hz 
QPO of SGR 1806-20 as U2, and the 92 and 150 Hz ones as the 
second and fourth overtones, one could associate one of the lower 
frequencies, 18 or 26 Hz, with Ul. 

• We have performed numerical simulations to test whether pure 
crustal shear oscillations can survive in the gaps of the Alfven con- 
tinuum, as it was proposed by van Hoven & Levin (2011). For 
this purpose we constructed a particular model (F) with a very flat 
spectrum which has very large gaps in between successive over- 
tones. In none of our simulations we could find the presence of 
crustal shear modes. Instead, predominantly Alfven oscillations of 
the continuum were present. This seem to contradict the results of 
Colaiuda & Kokkotas (201 lb), van Hoven & Levin (201 1) and van 
Hoven & Levin (2012). The mismatch with the results of Colaiuda 
& Kokkotas (2011b) may be due to a different interpretation of 
simulation results. In the work of van Hoven & Levin (2011) and 
van Hoven & Levin (2012) the assumption is made that the oscil- 
lations in the crust form a discrete spectrum, which coincides with 
the non-magnetized crustal spectrum (van Hoven & Levin 2011) 
or a magneto-elastic discrete spectrum (van Hoven & Levin 2012). 
This ansatz is a very good approximation for weak magnetic fields, 
when the crust still behaves approximately like a normal mode sys- 
tem. However, we think this assumption does not hold, in general, 
in the presence of a strong magnetic field, where the oscillations in 
the crust can no longer be described by an eigenvalue problem. 

• For mixed poloidal-toroidal magnetic field configurations 
(type I) where the toroidal field is confined in the region of closed 
poloidal magnetic field lines inside the star, we estimate the part 
of the Alfven spectrum of axial oscillations in the region of open 
field lines, where the toroidal component vanishes. Qualitatively, 
the spectrum is similar to the one of purely poloidal fields (type II) 
in that region. With increasing field strength of the toroidal compo- 
nent, the poloidal field strength near the closed field lines increases, 
too. This also means that the surface field becomes stronger. Ta- 
ble 6 (3rd set) shows the block for the different configurations bo 
we see that the necessary equivalent dipole magnetic field strength 
to match the observed QPO frequencies in the spectrum increases 
with increasing toroidal field, because of the stronger increase of 



9 30, 92, and 150 Hz for SGR 1806-20, and 28, 84, and 155 Hz for SGR 
1900+14 



the poloidal field strength close to the surface, while the field in the 
center is only mildly effected. 

• Configurations with an additional quadrupole-like component 
have a richer spectrum when this component is of comparable 
strength as the dipole-like one. For Q/D = 0.5, we find three 
turning points in the spectrum. The additional turning points can 
produce oscillations which could be used to match a larger number 
of observed frequencies. For a quadrupole dominated configura- 
tion Q/D > 1.0, we find asymmetric oscillations in both hemi- 
spheres of the neutron star. These provide a source for additional 
QPO frequencies which could be matched to observations. Mod- 
els with Q/D >> 1, i.e. dominantly quadrupole-like configura- 
tions, are currently not favoured, since they require magnetic field 
strengths lower than standard spin-down estimates. The presence of 
a large quadrupole-hke component in magnetars, comparable to the 
poloidal component inside the star is not ruled out by our model. 
In fact the possible existence of this strong component could be re- 
lated to the amplification of magnetic fields due to dynamo action, 
in initially rapidly rotating neutron stars in the standard magnetar 
formation scenario (Duncan & Thompson 1992). Once the dynamo 
stops, the magnetic field evolves mainly due to Hall drift (Goldre- 
ich & Reisenegger 1992) on a characteristic timescale 



THall = 5 10' 



Vlkmy 



(35) 



where B is the typical magnetic field strength in the magnetar in- 
terior, and L is the typical length scale of the magnetic field loops. 
Under such conditions, a quadrupole-like component of similar 
strength as the dipole-like one could still be present in magnetars, 
if it was present at birth. 

• Finally, we investigated magnetic field configurations limited 
to the crustal region. These configurations could be realized if the 
neutron star core is a superconductor of type I. In this case the shear 
modes cannot be damped into the core of the neutron star and they 
might survive sufficiently long to become observable. However, our 
simulation results suggest that these models provide no viable ex- 
planation for the observed frequencies in SGRs. We showed that 
the influence of the magnetic field strength on the frequencies of 
the magneto-elastic QPOs can be cast into a semi-empiric formula, 
which also contains the effect of different shear moduli. With this 
relation, we are able to scan a wide range of parameters of magnetic 
field strength and magnitude of the shear modulus. The two possi- 
ble types of QPOs are not favoured because: i) On one hand, the 
structure and in particular the frequencies of predominantly shear 
modes in the zero-magnetic-field case change significantly with in- 
creasing magnetic field, such that the lowest observed frequencies 
of / ^ SOHz cannot be reached for the currently expected val- 
ues of the shear modulus and the magnetic field strength, ii) On 
the other hand, predominantly Alfven oscillations will be quickly 
damped by phase mixing due and strong gradients in the spectra 
and an absence of turning points. For different EoS we obtain sim- 
ilar results. As expected, the core EoS does not play a large role in 
the determination of the critical magnetic field strength at which the 
frequencies begin to increase significantly with the magnetic field. 
The NV crust EoS leads to larger crusts with stronger shear moduli 
compared to the DH EoS. This, in turn, causes the necessary mag- 
netic field to change the predominantly shear QPOs significantly to 
be up to 2 X 10^^ G stronger than for the DH EoS. For quadrupole- 
or octupole-Uke dominated fields, the increase in frequencies of 
predominantly shear QPOs is less strong compared to dipole-like 
fields. However, from the neutron star spin-down measurement one 
has estimates of the dipole-like component of B15 ~ 1. At this 
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dipole-like field strength our previous arguments already apply and 
the frequencies are shifted to / > 30 Hz. 

To summarize our findings, we have explored a large parame- 
ter space of possible magnetic field configurations of neutron stars. 
We discovered additional features, such as new turning points in the 
spectra for purely poloidal configurations, which have the potential 
to explain more of the observed QPO frequencies. Configurations 
with a toroidal component, which is confined to the closed poloidal 
field lines, have very similar spectra compared to those of purely 
poloidal fields. Furthermore, our results do not favour the assump- 
tion that the magnetic fields are confined to the crust. 

Our magnetic field configurations do not include fields which 
have a toroidal component extending throughout the whole neu- 
tron star, which may be realized during a giant flare where strong 
toroidal fields are expected in the exterior (Beloborodov 2009; 
Fernandez & Thompson 2007; Nobili et al. 2008). For such con- 
figurations Colaiuda & Kokkotas (2011a) find discrete Alfven os- 
cillations. However, these fields require currents at the surface of 
the neutron star, which are difficult to model (and are set ad hoc in 
other studies). We postpone such considerations to future work. 

In a previous work (Gabler et al. 2012) we have studied the ef- 
fect of different EoS for one particular magnetic field configuration. 
The differences in the estimated magnetic field strength required to 
match observations that are caused by changing the EoS are of the 
same order (factor of a few) as the differences caused by assum- 
ing different magnetic field configurations. One of the main open 
questions in order to match observations is the effect of superflu- 
idity in the core of the neutron star, which probably leads to lower 
(and thus probably more realistic) estimates of the magnetic field 
strength. Furthermore, a model of how the coupling of the interior 
oscillations to the magnetosphere can lead to a modulation of the 
X-ray emission needs to be developed. We plan to address these 
two points in future studies. 
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APPENDIX A: MAGNETIC FIELD CONFIGURATIONS 
CONFINED TO THE CRUST 

A more complete derivation of the following can be found in Aguil- 
era et al. (2008). For axisymmetric fields it is possible to make the 
following ansatz 



where ^ = ^(r, 9) is a scalar whose angular part can be expanded 
according to 



^(r,^) = C^^^^ft(r), 



(A2) 



with C being a normalization constant. The poloidal magnetic field 
can thus be expressed as 



Br 



Be 



(A4) 



where x = r /rg and 5 is a constant related to C. For dipolar fields 
C = rsB/2. The (/^-component of the force-free condition V x 
B = ii'B leads to a Riccati-Bessel equation for Si (x) (see Aguilera 
et al. 2008). /i is a parameter related to the currents maintaining the 
magnetic field. For different /, we have 



Si{x) = aiiirsxji{iirsx) + bi/j.rsX ni{/j.rsx) , 



(A5) 



where ai and bi are parameters, and ji and ni are the spherical 
Bessel functions of first and second kind. The parameters ai and hi 
can be obtained by matching to an exterior solution of the magnetic 
field. 

For a given / both have to be determined by the boundary con- 
ditions for Br and Be at the surface. For a dipole external field we 
get 



ai — cos iiTs , 
bi = sin jjLTs , 

while for the quadrupole external field 



(A6) 
(A7) 



a2 



3^ ^ - sin//rs - 2— — + 3^ , (A8) 



^ COS fir s , ^sm//rs , ^smfiTs 

b2 = -3 , +cos/xrs-2 , ^, +3 , . (A9) 
(/irs)^ [firs) [firsy 

The octupole field can be matched by 

sin/xrs cos/xrs , sin/xrs 
as — cos firs — o lo- ^+30- 



firs 



(/irs)2 ifirs)^ 



+30 



cos firs 



. ^cosfirs sm/irs cos firs 
03 = smfMrs-\-5— ^ 15- ^—30- 



(firs) {firsY 



+30 



sm firs 



(AlO) 



(All) 



It remains to determine the value of fi. For magnetic fields 
confined to the crust = at the crust-core interface. This trans- 
lates into the condition Si — 0' ^nd hence for Z = 1 into 



tan [fi(rc 



fir cc - , 



(A12) 



Bpoioidai = V X (r X VV^) , 



(Al) 



which has to be solved numerically. For the quadrupole- and 
octupole-like configurations the corresponding equations are more 
complicated and are not given here. 

The particular form of A{x) in general does not have to co- 
incide with the solution in terms of spherical Bessel functions pre- 
sented above. But because of arbitrariness it is advantageous to use 
an analytical description here. 
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